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ABSTRACT 

The lensing of CMB photons by intervening large-scale structure leaves a characteristic imprint on 
its arcminute-scale anisotropy that can be used to map the dark matter distribution in projection on 
degree scales or ~ 100/i~^ Mpc comoving. We introduce a new algorithm for mass reconstruction which 
optimally utilizes information from the weak lensing of CMB anisotropics in the damping tail. Individual 
degree-scale mass structures can be recovered with high signal-to-noise from a foreground-free CMB map 
of arcminute scale resolution, specifically with a FWHM beam < 5' and a noise level < 15 (10~^-arcmin) 
or 41 (/iK-arcmin). 

Subject headings: cosmic microwave background - dark matter — large scale structure of universe 



1. INTRODUCTION 



It is well-known from the study of the weak gravita- 
tional lensing of faint galaxies that the distortion of back- 
ground images can be used to map the interve ning mass 



Squire^ 1993). 
Die, 



distribution in projection (Tyson et al. 1990; Kaiser & 



As the most distant background image 
available, maps of the CMB temperature distribution pro- 
vide a unique opportunity to map the distribution of dark 
matter. They provide information about structure s on the 
largest linear s cales in the high redshift universe ( |Zaldar- 
riaga & Seljak 1999) and hence complement information 
from weak lensing surveys. The main difficulty is that 
unlike an image of background galaxies, the temperature 
distribution of the CMB is to good approximation a Gaus- 
sian random field with no characteristic shape. 

Algorithms in the literature for extracting the interven- 
ing mass distribution from lensed CMB maps have shown 
the p otential for st atistical detections by the Planck satel- 
lite.!^ P^J^^Q-^deau] (1998) considere d the distortion to the 



Hessian of the temperature field. [Zaldarriaga fc Seljak 
(1999) considered distortions to the product of gradients 
of the temperature field. In neither case is it possible to 
extract hig h signal-to-no ise maps of the dark matter. 

Recently Zaldarriaga (2000) sho wed that the damping 
tail of CMB anisotropics (see e.g. |Hu fc Whit^ 1997) ex- 
hibits enhanced lensing effects in the four-point function. 
Indeed even the two-point function or power spectrum 
shows enhanced effects in this region due to the multi- 
tude of acoustic peaks and the sharp decline in intrinsic 
power associated with damping (Metcalf & Silk 1997). 

Hu (2001) showed that there is a quadratic estimator 
that recovers all of the information in the four-point func- 
tion about the mass distribution on large scales. Even 
for experiments like Planck that only partially resolve the 
damping tail, this estimator reduces the noise variance 
of the recovered projected mass power spectrum by over 
an order of magnitude. Planned experiments to measure 
arcminute-scale secondary CMB anisotropics can poten- 
tially use this statistic to map the dark matter at high 
signal-to-noise. 

We begin in by reviewing the effect of lensing on CMB 
temperature maps. In §^ we describe the reconstruc- 

^ http: / / astro. estec.esa. nl / S A-general/Projects/Planck 



tion algorithm and test it with realizations of the CMB 
temperature field and instrumental noise. We discuss ob- 
servational strategies for optimizing the reconstruction in 
For illustration purposes we use a fiat ACDM cos- 
mology throughout with parameters Vtc = 0.3, = 0.05, 
= 0.65, h = 0.65, n = 1 said Sh = 4.2 x 10"^. 

2. LENSING 

Weak lensing of the CMB photons by the intervening 
mass distribution remaps the primary temperature field 
0(n) as a func tion of the directional vec tor n on the sky 
as (e.g. ^eljakl 1996; [Goldberg fc Spergel| 1999) 



e(n) =e(n + V(/)), 



(1) 



where is the defiection angle which is related to the 
gravitational potential ^(x, I^) as 
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where D is the comoving coordinate distance along the line 
of sight and Da is the comoving angular diameter distance 
associated with D. Ds is the coordinate distance to the 
last scattering surface. In a fiat universe Da = D. The 
projected potential (/)(n) has a power spectrum Cff^ which 
is itself a projection of the gravitational potential power 
spectrum (see e.g. Hu 2000, Eqn. 28). For the scales of 
interest here, the gravitational potential power spectrum is 
in the linear regime and hence (j){n) is a Gaussian random 
field. The power spectrum of the defiection angles is given 
by L{L + 1)C^^ (see Fig. |l| for ACDM). For reference, the 
convergence power spectrum, or projected mass, is given 
by [L{L + l)]2C^^/4. The rms defiection angle 
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(3) 



is 2.6' for this model but its coherence scale corresponds 
to the peak of the log power spectrum at L ~ 60 or a 
few degrees. Counterintuitively then, the arcminute scale 
structure of the CMB temperature field yields information 
about the mass distribution on much larger linear scales. 



Fig. 1. — Lensing deflection power spectrum for the ACDM model. 
Error bars represent the total (sample plus noise) variance and sam- 
ple variance from recovery from an area of /sky =0.1 and an experi- 
ment with a beam of cr = 1.5^ and noise w~^l^ = 10 (10~^-arcmin). 
Note that the errors for L < 200 are dominated by sample variance 
implying that the recovered map has high signal-to-noise. 




The deflection power comes mainly from structures at red- 
shifts z 1 — 2 and scales of k ^ few xlO~^ Mpc~^ (see 



Zaldarriaga fc Seljakl 1999, Fig. 7). 

To simulate a lensed CMB map, one makes a Gaussian 
random realization of the unlensed CMB power spectrum 
Ci and remaps the temperature field according to a ran- 
dom realization of the projected potential Cf^^. Detector 
noise and residual foregrounds are then added as a real- 
ization of (7[^oise detector noise and a finite beam of 



a (FWHM) (Knox 1995) 
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-lg^(Z+l)crV81n2 



(4) 



with is the noise in units of (AT/T-radian)^. 

Because the defiection angles are small compared with 
the scale of the structures, the lensing effect is difficult 
to see directly in a map. To gain a better intuition for 
the nature of the effect, let us first consider lensing by a 
circularly symmetric Gaussian profile in projected mass 
with a scale of 5° and an amplitude corresponding to a 
3° maximum defiection. As in the case of the weak lens- 
ing of faint galaxies, the distortion represents a tangential 
shearing of the image. Unlike the faint galaxy case, the 
source image is a Gaussian random field. Although the 
temperature map itself clearly shows evidence for lensing, 
the two point statistics of the CMB temperature field do 
not suffice to reconstruct the mass distribution of the lens. 

3. RECONSTRUCTION 

The case of the symmetric lens in Fig. ^ suggests that 
a statistic related to the Laplacian of the temperature 
field would trace the underlying mass distribution. The 
fact that both hot and cold spots are l ense d alike washes 
out the signal in the Laplacian itself. Hu (2001) showed 
that a related statistic, the divergence of the temperature- 
weighted gradient of the map retains all of the inf orma- 
tion inherent in the four-point function ( ^aldarriaga 2000). 
Here we consider its use in mapping the dark matter. 




Fig. 2. Top: A 32° x 32° realization of the CMB temperature 
field. Bottom: A toy example of the lensing effect. A circularly 
symmetric projected mass with defiection angles comparable to the 
size of the structure. The distortion of the fine-scale anisotropy of 
the CMB traces the lensing structure on much larger scales. 



We describe the technique for reconstructing the dark 
matter field on small sections of the sky ^map 

< 60° 

for which spherical harmonic analysis can be replaced by 
Fourier techniques. For the generalization to the curved 
sky see^ (2001). 

The first step is to take the gradient of the temperature 
map by filtering in the Fourier domain. 



(5) 



where C,*°* = Q + C,"°"^^ Q (Ci) is lensed (unlensed) 
power spectrum. Note that taking the gradient effectively 
high-pass filters the map. Next construct an explicitly 
high-pass filtered temperature map 



Win) = I 

to weight the gradient 
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G(n) = VF(n)G(n) . 



(6) 



(7) 



3 




so that 



Fig. 3. Top: A 32° x 32° realization of the deflection fleld in 
the ACDM model. Bottom: Recovery of the deflection fleld with a 
(7 = 1.5' beam (FWHM) and detector noise of w'^/'^ = 10 (IQ-^- 
arcmin). 



Finally take a filtered divergence of this field in the Fourier 
domain 



D{h) 



J (27r)2 L 



iL ■ G(L)e*^-" . 



(8) 



The normalization factor N^/L may be chosen so that 
D{ri) averaged over an ensemble of CMB realizations re- 
covers the defiection field 



d(n) 



(27r)2 



(9) 



To determine Nl consider the operations directly in the 
Fourier domain, 



Nl f dHi 



(L.lia+L.l2a)jAlg|_, (10) 

where I2 = L — li. Ta ylor expanding Eqn. (|l|) for the 
lensing one obtains (Hu 2000) 



6(1) = 6(1) - 1 ^eihmi - ii)(i - ii) • ii , (11) 
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Notice that the filters are designed so that the lensing ef- 
fects in the Fourier domain add coherently such that the 
dot product above comes in as the square. This is a re- 
fiection of the optimization. 

We show an example of this reconstruction on a 32° x 32° 
field in Fig. (||) with detector noise added as appropriate 
for a beam of a = 1.5' and noise of w~^^'^ = 10 (10~^- 
arcmin) or (27/iK-arcmin) additionally low pass filtered 
to show L < 150 where the signal-to-noise is the highest. 
Alternately, Weiner filtering can be used to get a better 
visual impression of the fidelity of the map. In any case, 
the degree scale features in the map are recovered at good 
signal-to-noise. 

In the idealization of an ensemble of CMB maps lensed 
by the same structure, D^L) returns an unbiased estimate 
of the deflection map. However given that we only have 
one realization of the lensing per lens, it is important to 
understand the properties of the noise introduced by the 
Gaussian primary anisotropics themselves and the instru- 
mental and/or foreground noise. Following |Hu| (2001), 

{D''{L)D{V)) = (2^)2(5(L - V) (^L^Cf + Nl) , (14) 

so that Nl also plays the role of the noise power spec- 
trum. A deflection power spectrum extracted from this 
statistic must remove this noise bias. As discussed in Hu 
(2001), the noise bias may alternately be eliminated by 
cross-correlating maps reconstructed from independent l- 
bands in the original lensed map. 

Under the assumption of Gaussian statistics, the signal- 
to-noise per L in the deflection power spectrum is given 

by 



S_y _ 2L + 1 



2 yL'^Ct'' + Nl^ 



(15) 



and the precision with which the binned deflection power 
spectrum can be recovered by an experiment with a sky 
fraction of /sky = 0.1 (- 4000deg2) and a = 1.5', ^^/^ = 
10 (10~^-arcmin) is shown in Fig. [l| Also shown are the 
errors provided by sample variance alone. Since sampling 
errors dominate at L <^ 150, the recovered map has good 
signal-to-noise on those characteristic structures. This is 
independent of the actual sky fraction covered by the ex- 
periment. 

4. DISCUSSION 

The statistic introduced here utilizes CMB structures 
in the arcminute regime of the damping tail to map the 
dark matter on degree scales. As with the weak lensing of 
faint galaxies, image distortions manifest on small angular 
scales are used to reconstruct the mass on a much larger 
scale. Mapping the dark matter distribution therefore re- 
quires high resolution, high signal-to-noise maps of the 
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- noise 
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1 10 
detector noise w-^'^ (lO-^-arcmin) 

Fig. 4. — Signal-to-noise as a function of detector noise, beam 
and sky fraction /sky Solid lines include both sample and Gaus- 
sian noise (primary CMB and instrumental) variance; dashed lines 
include only Gaussian noise variance. The signal-to-noise drops off 
rapidly for w~^^'^ > 15 (10~^-arcmin) and FWHM beams a > 5^ 



CMB anisotropies themselves. Conversely, though a wide 
field of at least several degrees on the side is required to 
map the full extent of the structures expected, the statis- 
tic essentially high pass filters the input CMB maps. A 
true map that retains correlations across these scales is 
not necessary. 

To see how an observing strategy might be optimized 
for mapping the dark matter, let us consider the trade- 
offs between sky coverage, instrumental noise and beam. 
Because this statistic is a quadratic function of the temper- 
ature fluctuation data, the balance differs from the usual 
case. In Fig. ^ we show the total signal-to-noise in the 
measurement of the deflection power spectrum (summed 
in quadrature over L) of an experiment as a function of 
these parameters. We consider separately the case of noise 
variance from the Gaussian random primary anisotropies 
and detector noise alone and combined with the sample 
variance of the lensing fields. When the former exceeds the 
latter, a high signal-to-noise map of the structures results. 
Because this is an integrated statistic, the characteristic 
signal-to-noise for large-scale features is much higher (see 
Fig. 0). 

Compare the steep increase in the signal-to-noise as the 
detector noise is reduced with the shallow increase with 
sky coverage of /gYy • Up until w~'^^'^ ~ 10 (10~^-arcmin), 
observing time is best spent going deep rather than wide. 
Beyond this point, the intrinsic noise variance provided by 
the primary CMB anisotropies themselves begins to domi- 
nate and saturate the signal-to-noise. If the goal is to pro- 
duce a high signal-to-noise map of structures, then going 



down to 



-1/2 



1 (10 -arcmin) can achieve substan- 



tially improved maps of the finer scale structures in the 
map. Another crucial factor is the beam size. To resolve 
the structures that best trace the lensing, a beam of a < 5' 
is required and it is not until a ~ V — 2' that the gains 
saturate. If foregrounds are not removed from the map 
through their spatial coherence and/or frequency depen- 
dence then this balance can shift to larger angular scales 
and more sky coverage. For the 1.5^ 10 (10~^-arcmin) 



baseline experiment, inclusion of Gaussian random noise 
from the Sunyaev-Zel'dovich and Vishniac effects in 
imply a relative degradation in signal-to-noise of ~ 10% 
and ~ 1% (for r = 0.1) respectively and so do not require 
substantial reoptimization. 

A high signal-to-noise map of the dark matter in pro- 
jection can also be used to pull out tracers of the large- 
scale structure of the universe in other maps through cross- 
correlation. Examples include secondary anisotropies such 
as th e intefflated Sachs-V\ ^olfe and Sunyaev-ZePdovi ch ef- 
fects (iGoldberg fc Spergel l999; |Seljak fc Zaldarriaga| l998; 
Cooray Hu| 2000). One can show that the statistic em- 



ployed here retains all of the information in the full bis- 
pectrum of the secondary-lensing-primary correlation and 
so is the optimal statistic to measure these correlations. 

The filters used to reconstruct the dark matter map for- 
mally require as input the power spectrum of the CMB 
before lensing. The lensed CMB power spectrum will of 
course be measured to exquisite precision by CMB satel- 
lites and by the input temperature map themselves. Em- 
ploying the lensed CMB power spectrum in the filter or 
an otherwise slightly incorrect assumption simply degrades 
the signal-to-noise by a correspondingly small amount but 
does not introduce spurious structures in the ensemble- 
averaged recovery. They appear as a calibration error for 
the mass map. Indeed in the context of a parameter- 
ized cosmology the unlensed CMB power spectrum may 
itself be reconstructed from the observed spectrum. For 
a non-uniform survey geometry, with perhaps foreground- 
contaminated regions removed, more sophisticated tech- 
niques than the Fourier-transform filtering scheme em- 
ployed here will have to be developed. These complications 
should not present an insurmountable obstacle to the goal 
of mapping the dark matter in projection at intermediate 
redshifts. 

Acknowledgements: I acknowledge useful conversations 
with A.R. Cooray and M. Zaldarriaga as well as support 
from NASA NAG5-10840, DOE OJI and an Alfred P. 
Sloan Foundation Fellowship. 



REFERENCES 

Bernardeau, F. 1998, A&A, 338, 375. 
Cooray, A.R. Hu, W. 2000, ApJ, 534, 533. 

Goldberg, D.M. Spergel, D.N. 1999, Phys. Rev. D, 59, 103002. 

Hu, W. 2000, Phys. Rev. D, 62, 043007. 

Hu, W. 2001, Phys. Rev. D, in press, astro-ph/0105117. 

Hu, W. White, M. 1997, ApJ, 479, 568. 

Kaiser, N. Squires, G. 1993, ApJ, 404, 441. 

Knox, L. astro-ph/9504054, Phys. Rev. D, 52, 4307. 

Metcalf, R.B. Silk, J. 1997, ApJ, 489, 1. 

Seljak, U. 1996, ApJ, 463, 1. 

Seljak, U. Zaldarriaga, M. 1999, Phys. Rev. D, 60, 043504. 
Tyson, J.A. Wenk, R.A. Valdes, F. 1990, ApJ, 349, 1. 
Zaldarriaga, M. 2000, Phys. Rev. D, 62, 063510. 
Zaldarriaga, M. Seljak, U. 1999, Phys. Rev. D, 59, 123507. 



